Quantum Ghost Imaging through Turbulence 



O 

(N 



Oh. 



> 
00 

m 
m 

O 



X 



p. Ben Dixon/ Gregory Howland/ Kani Wai Clifford Chan,^ Colin O'SuUivan-Hale,'^ Brandon Rodenburg,'^ 

Nicholas D. Hardy,"' Jeffrey H. Shapiro,'' D. S. Simon,^ A. V. Sergienko,^ R. W. Boyd,^'^ and John C. Howeh^ 

'Department of Physics and Astronomy, University of Rochester, Rochester, New York 14-627, USA 
'^Rochester Optical Manufacturing Company, 1260 Lyell Avenue, Rochester, NY 14-606, USA 
'^Institute of Optics, University of Rochester, Rochester, New York, 14627, USA 
"^Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 
^Department of Electrical & Computer Engineering, Boston University, Boston, MA 02215, USA 

(Dated: January 12, 2013) 

We investigate the effect of turbulence on quantum ghost imaging. We use entangled photons 
and demonstrate that for a novel experimental configuration the effect of turbulence can be greatly 
diminished. By decoupling the entangled photon source from the ghost imaging central image 
plane, we are able to dramatically increase the ghost image quality. When imaging a test pattern 
through turbulence, this method increased the imaged pattern visibility from V = 0.14 ± 0.04 to 
V = 0.29 ± 0.04. 

PACS numbers: 42.68.Bz, 42.30.Va, 03.67.Hk 



Introduction. — The phenomenon of ghost imaging 
(GI), first observed by Pittman et al. in 1995 1], is 
a method of generating the image of an object from 
correlation measurements. Pittman's experiment made 
use of pairs of entangled photons. One of the photons 
passed through a transmission object and then to a pho- 
ton counter with no spatial resolution. The other photon 
passed directly to a spatially resolving photon counter. 
When looking at coincident photon detections, the de- 
tectors were able to see the object despite the fact that 
the object and the spatially resolving detector were in 
different arms of the experiment. While it was initially 
thought to be a quantum mechanical effect reliant upon 
the entanglement between the two photons, similar re- 
sults were later obtained using classical sources 0]. 

In addition to clarifyirig the boundary between quan- 
tum and classical effects [3|-l5| , GI has been used for lens- 
less imaging ^ , super-resolution imaging 0, Q , and en- 
tanglement detection [9] . More recently, research has rec- 
ognized connections between GI and compressive sensing 
llj . The distributed nature of GI has made it a 



Spatially- 
resolving 
Detector 




Image Arm 



NLC 



Correlator 
1 O • I 



Lens 



50/50 
BS 



Bucket 
Detector 



Laser Filters Lens 
Object Arm — - 



Object; SPAD i 



[1 

candidate for distributed image processing, and for dis- 
tributed sensing and communication schemes. 

For many optical applications, imaging through tur- 
bulence is unavoidable 
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GI is no different and 
the effect of turbulence on GI performance has begun to 
be theoretically investigated 12l4l4|. In this paper, we 
experimentally investigate the effect of turbulence on GI 
using entangled photons. We introduce a novel GI setup 
that allows us to minimize the effect of turbulence. 

Theoretical description. — The experimental apparatus 
is depicted in Fig. [TJA biphoton state IV') is created at 
a nonlinear crystal |l5| and then split by a 50/50 beam- 
splitter, sending the biphoton into two arms of the appa- 
ratus. 

In the object arm, the biphoton travels a distance 
2/ + A, to a lens which has focal length /. The bipho- 



FIG. 1: (color online). A 325 nm pump beam undergoes 
SPDC at a nonlinear crystal (NLC), the output passes a 
beamsplitter (BS). One beam is sent through a lens and onto 
a transmission object. The other beam is sent through a lens 
and onto a scanning slit. The ghost image of the object is 
profiled by the slit. Photons are detected with single-photon 
avalanche diodes (SPAD). A coincidence circuit correlates the 
measurements. 



ton then travels a distance 2/ to a photon detector with 
no spatial resolution (a "bucket" detector). A transmis- 
sion object — consisting of alternating opaque and clear 
vertical bars — is placed just in front of the detector. 

In the image arm, the biphoton travels a distance 2/ — 
A to a lens which again has focal length /. The biphoton 
then travels a distance 2/ to a spatially-resolving photon 
detector. 

For A = the detectors and crystal are all located at 
image planes of each other. As one arm's lens/detector 
is moved towards the crystal by a distance A, the other 
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FIG. 2: (color online). The experimental configuration is 
shown conceptually using the Klyshko picture [l6[, the object 
mask (on the right) is ghost imaged onto the scanning slit 
(on the left). The nonlinear crystal is offset from the cen- 
tral image plane by a distance A. The top picture shows the 
turbulence — represented by wavy lines — between the crystal 
and the lens. The bottom picture shows the turbulence lo- 
cated between the lens and the object. Experimentally rele- 
vant distances are labelled. 



arm's lens/detector is moved away by the same distance, 
keeping the sum of the arm's length constant, see Fig. [2] 

Turbulent air flow is introduced into a thin region in 
the beam path of the object arm. For turbulence be- 
tween the crystal and the lens, it is a distance from 
the crystal — or a distance — A from the central image 
plane. For turbulence between the lens and the object, 
it is a distance A — Zi from the object. 

The relevant function for GI is the second order de- 
gree of coherence G^'^\xi,X2), where xi is a transverse 
position variable in the plane of the spatially-resolving 
detector and X2 is a transverse position variable in the 
plane of the bucket detector. We begin with the standard 
quantum mechanical form and include an additional en- 
semble averaging — represented by large outer brackets — 
to account for the statistical effect of turbulence: 



G'^^\xi,X2) = ({^\^{xl)El{x2)Es(x2)%{xlm 



(1) 

Neglecting overall normalization, this can be represented 
in the following way: 
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(2) 

Subscript s and i indicate variables in the crystal plane 
and subscript t indicates variables in the plane of the tur- 
bulence. The function ip(xs,Xi) is the transverse bipho- 



ton wavefunction which we approximate as a plane- wave 
with delta function correlations 'il){xs,xi) = 6{xs — Xi). 
The function 'K{xs,X2]Xt) is a propagation operator go- 
ing from the crystal plane to the object arm detection 
plane, passing through the plane of turbulence; H(xi,xi) 
is a propagation operator going from the crystal plane to 
the image arm detection plane. These operators can be 
represented in the following way: 
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We are assuming a narrow sheet of turbulent air, whose 
effect on propagation can be characterized by a multi- 
plicative operator T{xt). In our theoretical treatment, 
we assume that the lenses are sufficiently large that they 
capture all of the light from the SPDC source. As a re- 
sult, both turbulence locations in Fig. [2] are governed by 
the same operators. 

We model the turbulence as a 6/3 scaling law effect: 
(T* {it)T{xt)') — exp [—a {xt — it)^/2] , where a param- 
eterizes the strength of the turbulence and has units 1 / m^ 
[17I ITsI l . After integration, the resulting expression for 
G(2)(a:i,X2) is: 
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The ghost image X(a:i) is then the product of the object 
and G(^)(a;i, 2:2), integrated over X2- We represent the 
object as: 0(2:2) = exp [— x|/2w^] (l-t-cos(A;oa:2))- Here 
w is the spatial width of the illuminating beam and ko is 
wavenumber for the object's pattern spacing. 

Assuming {li — A)y^ <^ kw, the ghost image is found 
to be 



2{xi) — exp 



li^T] (l + Vcos(fc„.r)). 



(6) 



I{xi) has the same form as 0{xi) with the object's 
unity visibility replaced by the detected ghost image vis- 
ibility V: 



V = g X exp 



a {h - A)^ 
2{k/kof 



(7) 



Where g is the the optimum ghost image visibility for 
turbulence. As the turbulence increases in strength (in- 
creasing a), the detected visibility V decreases — thus ob- 
scuring the detected pattern. Similarly, when the turbu- 
lence is moved away from either the central image plane 
or the detector, li — A increases, thereby decreasing V. 
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Experiment. — Collimated light from a 3 mW, 325 nm 
HeCd laser with a 1/e^ full width of approximately 
1600 fim pumped a 10 mm thick BBO nonlinear op- 
tical crystal. The crystal was oriented for degenerate 
type-I collinear spontaneous parametric down-conversion 
(SPDC). After the crystal, the pump beam was blocked 
by colored glass filters and the SPDC bandwidth was 
limited by a 3nm wide spectral filter centered at 650 nm. 
The remaining SPDC beam was split into two arms by a 
50/50 beamsplitter. 

In the image arm, a lens was located 1000 mm — A 
from the crystal; in the object arm, a lens was located 
1000 mm + A from the crystal. Both lenses had focal 
length / — 500 mm. Detectors were located 1000 mm 
from the lenses. 

The transmission object was a test pattern located 
1000 mm from the lens. The bucket detector consisted of 
a lOx microscope objective which collected the transmit- 
ted light into a multimode optical fiber. The pattern had 
unity visibility and 3.6 cycles per mm, which resulted in 
an object pattern wavenumber of ko = 7.2 x tt mm^^ 

The spatially-resolving detector consisted of a com- 
puter controlled scanning slit located 1000 mm from the 
lens, which was again followed by a 10 x microscope ob- 
jective which collected the transmitted light into a mul- 
timode optical fiber. The slit had an approximate width 
of 40 /xm and was scanned in 5 /zm increments, giving 
spatial resolution. 

The optical fibers were connected to Perkin Elmer 
single-photon avalanche diode detectors. The outputs 
of these detectors were time correlated using a PicoHarp 
300 from PicoQuant. Photon counts were integrated at 
each slit location for between 1 and 4 seconds. The spa- 
tially resolved coincident detections made up the ghost 
image profiles. 

A downward pointing heat gun was mounted above the 
setup, providing turbulent air flow across the beam path. 
The effect of the turbulence was fitted to the model's 
wave structure function a [3| ■ From the fit we deter- 
mined a = 2.5 ± 1.5 mm~^. 

Data was taken for an unshifted configuration with 
A = 0, and for a shifted configuration with A = 330 
mm. In each configuration, ghost images were recorded 
with turbulence present in the object arm: both between 
the crystal and lens, and between the lens and the object. 
Ghost images were also recorded with no turbulence. 

While allowing access to the central image plane of the 
apparatus, the shifted configuration introduced two ex- 
perimental limitations: the flux through the object and 
scanning slit decreased significantly as a result of the de- 
tectors being away from the beam focus, and fewer spatial 
frequencies contributed to the ghost image as a result of 
the nonlinear crystal having a stronger aperturing effect. 

The recorded ghost image profiles were fitted to T{xi) 
from Eq. [Bl The fit included a visibility term which con- 
stituted our measurement of the visibility V. 
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FIG. 3: (color online). Representative ghost images for 
the unshifted configuration (left), and shifted configuration 
(right). The top row shows images with no turbulence. The 
middle row shows images for turbulence between the lens and 
the object, 203 mm (right) and 229 mm (left) from the ob- 
ject. The bottom row shows images for turbulence between 
the crystal and the lens, 432 mm from the crystal. Points are 
experimental data while cuves are fits to the data. Counts 
are measured in coincident photon detections per second. 



Representative ghost images are shown in Fig. |31 
With no turbulence, the unshifted configuration pro- 
duced ghost images, with visibilities of 1.00 ±0.05. With 
no turbulence, the shifted configuration produced ghost 
image visibilities of only 0.65 ±0.05. The scans also show 
the decreased flux and the broader beam profile associ- 
ated with the shifted configuration. 

Visibilities for turbulence between the lens and the ob- 
ject are shown in Fig. |31 For both the shifted and un- 
shifted configurations, when the turbulence was close to 
the object, the observed visibility was near its no turbu- 
lence level. As the turbulence was moved away from the 
object, the ghost image visibility decreased, thereby de- 
grading the image. The visibility for the unshifted con- 
figuration remained above the visibility for the shifted 
configuration for all turbulence locations. 

Visibilities for turbulence between the crystal and the 
lens are shown in Fig. [S] This is the main result of the 
experiment. Visibilities decreased as the turbulence was 
moved away from the crystal, however, the unshifted con- 
figuration had lower fringe visibility than the shifted con- 
figuration. Indeed, for turbulence located 482 mm from 
the crystal, the visibility was V = 0.14 ± 0.04 for the un- 
shifted configuration, while for the shifted configuration 
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FIG. 4: (color online). Ghost image visibilities are shown for 
turbulence between the lens and the object. Visibilities are 
plotted as a function of distance from the object to the turbu- 
lence {h — A in Eq. [Tjl. Data for the unshifted configuration 
are shown as blue circles while data for the shifted configu- 
ration are shown as purple squares. Curves are plots from 
Eq. [T] The solid curve is for the unshifted configuration, with 
g = 1.00, and A = 0. The dashed curve is for the shifted con- 
figuration, with g = 0.65 and A = 330 mm. For both curves 
a = 2.0 mm"^. 
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FIG. 5: (color online). Ghost image visibilities are shown for 
turbulence between the crystal and the lens. Visibilities are 
plotted as a function of distance from the crystal to the tur- 
bulence (li in Eq. [T]). Data for the unshifted configuration 
are shown as blue circles while data for the shifted configu- 
ration are shown as purple squares. Curves are plots from 
Eq.[7l The solid curve is for the unshifted configuration, with 
g = 1.00, and A = 0. The dashed curve is for the shifted 
configuration, with g — 0.65 and A — 330 mm. For both 
curves a = 2.0 mm~^. The vertical line marks the location of 
the central image plane. 



it was V = 0.29 ± 0.04. Moving to the shifted configura- 
tion doubled the visibility. Data was taken for turbulence 
as close as 380 mm to the crystal, or 50 mm from the cen- 
tral image plane. 

Concluding remarks. — By moving the crystal from the 
central image plane we were able to place turbulence in 
this plane. This decreased the observed effect of turbu- 
lence, in fact it more than made up for the inherent loss of 



visibility associated with the shifted configuration. This 
technique has use in free space GI applications where 
turbulence is involved. By arranging detectors to place 
an image plane at the location of the turbulence, image 
degradation from the turbulence can be diminished. 

Although we have used entangled photons, similar re- 
sults are expected for thermal light GI. It should also 
be noted that the theoretical description assumes delta 
function correlations for the biphoton state and uses a 
thin region, non-Kolmogorov turbulence model 17- 19|. 
We are currently extending our theoretical description 
to include different SPDC correlation areas and a more 
complex description of turbulence including the possibil- 
ity of volume turbulence. This will be presented in a 
forthcoming paper. 

In this paper we have demonstrated a method of ame- 
liorating the effects of turbulence on GI systems, and 
have provided a theoretical model which accurately de- 
scribes the experimental data. We shift the source of en- 
tangled photons away from a quantum GI system's cen- 
tral image plane, and place turbulence near this plane. 
This dramatically increases the ghost image contrast. 
For turbulence located 482 mm from the crystal, this 
technique took the observed pattern visibility from V = 
0.14 ± 0.04 to V = 0.29 ± 0.04, doubling the system's 
imaging visibility. 

We acknowledge discussions with J. H. Eberly and sup- 
port by DARPA DSO InPho grant W911NF-10- 1-0404 . 
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